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Abstract 

. The generating function of double Hurwitz numbers is known to 

^P^l become a tau function of the Toda hierarchy. The associated Lax 

1^ ' and Orlov-Schulman operators turn out to satisfy a set of generalized 

■ string equations. These generalized string equations resemble those 
of c = 1 string theory except that the Orlov-Schulman operators are 
contained therein in an exponentiated form. These equations are de- 

CN ' rived from a set of intertwining relations for fermion bilinears in a 

^ . two-dimensional free fermion system. The intertwiner is constructed 

l/-^ [ from a fermionic counterpart of the cut-and-join operator. A classical 

■ limit of these generalized string equations is also obtained. The so 
called Lambert curve emerges in a specialization of its solution. This 

I seems to be another way to derive the spectral curve of the random 

matrix approach to Hurwitz numbers. 
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1 Introduction 

Hurwitz numbers count the topological types of finite ramified coverings of a 
given Riemann surface [Tj. Some ten years ago, Hurwitz numbers of coverings 
of the Riemann sphere CP^ turned out to be related to Hodge integrals on 
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the Deligne-Mumford moduli space Mg^n of marked stable curves [21 El S] 
and Gromov-Witten invariants of CP^ 0, El [H]- These observations led 
to various developments and applications, such as new proofs [HI |12l |13] 
of Witten's conjecture [9] (Kontsevich's theorem [TO]) on two-dimensional 
topological gravity, in other words, integrals of ■j/'-classes on Adg^n- Recently, 
a set of new recursion relations for Hurwitz numbers were derived [HI HSl [IS] 
as an analogue of Eynard and Orantin's "topological recursion relations" [T7] . 

These developments in the last decade are also more or less connected 
with integrable hierarchies of the KdV, KP and Toda type (see reviews by 
the Kyoto group [HI [191 120] for basic knowledge on those integrable systems) . 
Firstly, several Toda-like equations show up in relation with Gromov-Witten 
invariants of CP^. Secondly, the Witten conjecture is formulated in the 
language of the KdV hierarchy and its Virasoro symmetries. Subsequent 
studies [211 [221 123 on the Witten conjecture from the Hurwitz side revealed 
a connection with the KP hierarchy as well. 

A master integrable system in this sense is the (two-dimensional) Toda 
hierarchy f24j. As pointed out by Okounkov a generating function of 
"double Hurwitz numbers" gives a special solution of the Toda hierarchy. 
Being a tau function of the Toda hierarchy, this generating function has two 
sets of independent variables t = {ti,t2,- ■ ■) and i = {ii,i2, ■ ■ ■)■ By spe- 
cializing part of these variables to particular values, the generating function 
of "simple Hurwitz numbers" are recovered. The aforementioned Toda-like 
equations and KP hierarchy may be thought of as equations satisfied by those 
specialized tau functions. 

The goal of this paper is to derive the generalized string equations 

for this special solution of the Toda hierarchy, and to examine implica- 
tions thereof. Here /3 and Q are parameters prepared along with t and 
i. L, L and M, M are, respectively, the Lax and Orlov-Schulman opera- 
tors, all being one- dimensional difference operators. Compared with the 
cases of two-dimensional quantum gravity [25j and c = 1 string theory 
[261 123 [2H1 (23 ISDl [Slj, these generalized string equations have a notable 
new feature. Namely, they contain the exponentials e^^'^ and e^*^ of the 
Orlov-Schulman operators. This is in sharp contrast with the conventional 
string equations that are linear in M and M (except for the the case of a 
deformed version of c = 1 string theory [3D], in which the string equations 
are quadratic in M and M). 

The exponentials of M and M are related to the "quantum torus algebra" 
in a two-dimensional free fermion system. The same Lie algebraic structure 
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plays an important role in the melting crystal model of five-dimensional gauge 
theory as well [521 123] • We borrow some technical ideas developed therein. 
The generalized string equations are derived from algebraic relations (referred 
to as "intertwining relations") among fermion bilinears. A clue therein is 
a special fermion bilinear that corresponds to the so called "cut-and-join 
operator" |51J . 

The generalized string equations also have a classical limit of the form 

where £, A4, £, Ai are "long- wave" or "dispersionless" limits [20] of the Lax 
and Orlov-Schulman operators. Conceptually, the classical limit amounts to 
the genus- zero part of Hurwitz numbers. Similarity with c = 1 string theory 
becomes even more obvious in this limit (which amounts to the genus-0 part 
of string amplitudes). The generalized string equations of c = 1 string theory 
[2Sl EZl ESI 121] have a classical limit of the form 

C = CM, C-^ = C-^M. 

Thus, roughly speaking, the linear terms therein are now replaced by 

the exponential terms e^^,e^^. Thanks to this similarity, one can apply 
the method developed for solving the generalized string equations of c = 1 
string theory [29] to construct a solution of the foregoing equations in the 
form of power series of t^'s and t^s. This construction simplifies when either 
tfc's or tfc's are specialized to particular values for which the tau function 
reduces to the generating function of simple Hurwitz numbers. Remarkably, 
we encounter therein the equation 

X = ye^ 

of the "Lambert curve" that lies in the heart of the new recursion relations 
[m [15], [T6|. This should not be a coincidence. Presumably, the classical 
limit of the generalized string equations will be another way to derive the 
"spectral curve" in the random matrix approach [13 ESI ES] • 

This paper is organized as follows. Section 2 reviews the notion of Hurwitz 
numbers and their generating functions. Schur functions and their special 
values are used to interpret the generating functions as tau functions. Sec- 
tion 3 prepares technical tools from a two-dimensional free fermion system. 
We introduce relevant fermion bilinears, and recall Okounkov's result on a 
fermionic representation of the generating function of double Hurwitz num- 
bers [6]. Section 4 presents the generalized string equations. We start from 
intertwining relations of fermion bilinears, and translate those relations to 
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the generalized string equations. Section 5 formulates the classical limit of 
the generalized string equations. The classical limit is first derived by heuris- 
tic consideration, then justified by showing that the associated tau function 
has an /i-expansion. Section 6 is devoted to solving this classical limit of the 
generalized string equations. The solution is constructed in much the same 
way as that of the generalized string equations for c = 1 string theory. The 
Lambert curve shows up in a specialization of this solution. 

2 Generating functions of Hurwitz numbers 

In this section, we use various notions and formulae on partitions. Young 
diagrams and Schur functions. They are mostly borrowed from Macdonald's 
book [37]. 

2.1 Hurwitz numbers of Riemann sphere 

Let us consider finite ramified coverings of CP^. Two coverings tt : F — t- CP^ 
and tt' : r' — )■ CP^ are said to be topologically equivalent if vr and tt's are 
connected by a homeomorphism : F — F' as vr = tt' o 0. Let [vr] denote 
the equivalence class of the covering vr. Note that F can be a disconnected 
surface. 

Given a positive integer d and a set of distinct points Pi, ■ ■ ■ , of CP^, 
there are only a finite number of nonequivalent d-fold coverings that are un- 
ramified over all points other than Pi, ■ ■ ■ ,Pr- Those nonequivalent coverings 
can be further classified by the ramification data of sheets over Pi, ■ ■ ■ , Pr- 
These data are given by conjugacy classes Ci, ■ ■ ■ , C,. of the d-th symmetric 
group Sd- 

A conjugacy class C of Sd is determined by the cycle type 

IJ- = {fJ-i,- ■ ■ ,fJ'i), fJ'i > ■ ■ ■ > fJ'i, 1/^1 = yUi H ^f^i = d 

of a representative a E Sd of C, e.g., 

cr = (1, ■ ■ ■ , /ii)(/ii + 1, ■ ■ ■ , /ii + /i2) ■ ■ ■ (/ii H h fii-i, ■■■ ,d), 

where (ji, ■ ■ ■ ,jm) denotes the cyclic permutation sending ji — )■ j2 — ^ ■ ■ ■ — ^ 
jm — ^ ji- The cycle type thus becomes a partition of d, and has another 
expression 

fi= (l"^i2"^2---) 

with the numbers rrii of z-cycles. Let C(/i) denote the conjugacy class deter- 
mined by a partition of d. 



4 



A cyclic permutation (ji, ■ ■ ■ ,jm) of length m represents a cyclic covering 
of degree m realized., e.g., by the Riemann sheets of (z — d)^!^ above the 
point a. We use the cycle type /i = (yUi, ■ ■ ■ as ramification data above 
a point of CP^ to show that the sheets above that point locally look like a 
disjoint union of / cyclic coverings of degree /ii, /i2, ■ ■ ■ , 

The Hurwitz number Hd{Ci, ■ ■ ■ , Cr), also denoted by H(i{n^^\ ■ ■ ■ , /x*-^^) 
where /i'^'^^'s are the cycle types of C^'s, is defined to be the sum 

Hd{Ci, ■ • ■ , Cr) = 



|Aut( 



TT 



of the weights l/|Aut(7r)| over all equivalent classes [vr] of possibly discon- 
nected coverings vr that are ramified over the r points Pi, ■ ■ ■ ,Pr with the 
ramification data Ci, ■ ■ ■ ,(7^. Aut(7r) denotes the group of automorphisms 
of vr. As the notation suggests, Hd{Ci, ■ ■ ■ , Cr) does not depend on the po- 
sition of Pi, ■ ■ ■ ,Pr. 

The Hurwitz numbers can be determined by a genuinely group theoretical 
method (Burnside's theorem This leads to the beautiful formula 

\\\=d ^ ■ / fc=i 

where the sum is over all partitions A of d. Let us explain the notations used 
in this formula. 

Firsly, dim A denotes the dimension dim Vx of the irreducible representa- 
tion {px, Vx) of 5*^ determined by A. In other words, dim A is the number of 
all standard tableaux of shape A, and can be calculated by the so called hook 
length formula 

dimA = — — , (2) 

n M^'^') 

where h{i,j) denotes the length of the hook cornered at the cell {i,j) G A. 
Note that partitions are identified with the associated Young diagrams. As an 
immediate consequence, dim A turns out to be symmetric under the transpose 
(or conjugate) A t— t- *A of the partition, namely, 

dim*A = dimA. (3) 

Secondly, fxifJ') denotes the value fx{C{fi)) of the class function 

MO = MZl|q (4) 
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on Sd evaluated at C — C{fj,). xx denotes the character Try^ px- The cardi- 
nahty |C(//)| of C(/x) C can be written as 

\C{i,)\^-, z.^HrriiU^K (5) 
i>i 

The values of fx{C) for the simplest two cases C = (1*^), (l''~^2) can be 
explicitly calculated as 

Ml'') = l, fx{l'-'2) = ^ (6) 



where 



2^ 



This number has another useful expression of the form 

«:a = 2 J] (z-j), (7) 
(«j)eA 

which implies, e.g., the anti-symmetric property 

= (8) 



2.2 Simple Hurwitz numbers 

Let us review the notion of simple Hurwitz numbers. They are the Hurwitz 
numbers such that the types of all but one ramification points Pi, - • ■ ,Pr 
are restricted to l'^~^2. The exceptional ramification point Pr+i can have an 
arbitrary cycle type /i. The Hurwitz numbers of this type 




are called simple Hurwitz numbers. Introducing an extra (finite or infinite) 
set of variables x = {xi,X2, ■ ■ ■ ), one can construct a generating function of 
these numbers as 

r=0 d=0 \n\=d 
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where p^'s are the monomials 

Pn PlJ.lPl-12 

of the power sums 

Pk = ^x'l, /c = l,2,---. 

i>l 

We now substitute ([9]) in the definition of Z{x) and use the Frobenius formula 

j2 ^^^^p, = (10) 

\^,\=d 

to rewrite the sum over /j, into a Schur function times numerical factors. 
The numerical factors partly cancel with some other factors. The generating 
function thus reduces to 

Z{x) = J2^e^^^/^Q^^\s,{x). (11) 

xev ' 

Z{x) turns out to be a tau function of the KP hierarchy [181 [IH]- To 
see this, we change the variables from x to the standard time variables t = 
{ti,t2,- ■ ■) of the KP hierarchy via the power sums as 




and consider the Schur functions sx{x) as functions of t. It is convenient to 
use Zinn- Justin's notation s\[t\ for the latter [39j . 

Actually, sx[t] can be redefined directly. Let us recall the Jacobi-Trudi 
formula 

sx{x) = det{hx^^,+,{x))l^^,, (12) 

where Aj's are parts of A (which are assumed to be equal to for i > n, 
namely, A = (Ai, ■ ■ ■ , A„, 0, ■ ■ ■ )), and hm{xys are the complete symmetric 
functions defined by the generating function 

oo 

m=0 i>l 
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Let hm[t] denotes hm{x) considered as a function of t. /i^M's can be redefined 
by the generating function 

oo / oo 

m=0 \k=l 

Consequently, sx[t] can be expressed as 

sx[t] = detih,^^,^,[t])l^^,. (13) 
One can thus replace sx{x) by sx[t] to obtain the generating function 

Z[t]=J2^e^^^/^Q\^\sx[t]. (14) 



Xev 



Identifying Z[t] as a KP tau function requires some more consideration. 
Several methods are known in the literature [121 EB [221 [23] ■ One can explain 
this fact in the context of the Toda hierarchy as well. A clue is the formula 

dim A r T / N 

^^ = sa[1,0,0,---] (15) 

that can can be derived, e.g., from the Frobenius formula (fTOj) by letting 
Pi = 1 and pk = for k > 1. Consequently, Z[t\ can be rewritten as 

Z[t] = J2^^^'^'Q^'^^x[t]sx[l,0,---]. (16) 

Xev 

This function is a specialization of the generating function Z[t,i] of double 
Hurwitz numbers introduced below. Z[t,i] is a tau function of the Toda hi- 
erarchy (or, rather, the two-component KP hierarchy |18j. because the lattice 
coordinate of the Toda lattice is absent here). It is well known [2^ that any 
Toda (or two-component KP) tau function is also a tau function of the KP 
hierarchy with respect to one of the two sets of variables. This implies that 
Z\t\ is a KP tau function. 



2.3 Double Hurwitz numbers 

Let us choose yet another point Pq of CP^ of an arbitrary ramification type 
fi in addition to the r + 1 points in the case of simple Hurwitz numbers. The 
Hurwitz numbers of this type 

H,m, i'-^2,...,i^~% ^.) = y: (^X {yY fx{f^)fxm) (IT) 

r \x\=d ^ ■ ^ 
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are called double Hurwitz numbers. To construct a generating function of 
these numbers, we introduce a new set of variables x = {xi, X2,- ■ ■), their 
power sums 



pk = Yl 



and their monomials 



PX = PX1PX2 ■ ■ ■ 



along with the variables in the case of simple Hurwitz numbers. Again, with 
the aid of the Frobenius formula (llUp . the generating function 

r—O d=0 \^\ = \il\=d J. 

can be converted to a sum over all partitions as 

Z{x, x) = Y, e^^^l^Q\^\sx{x)sx{x). (18) 

AGP 

We now introduce the two sets t = (ti, ^2, ■ ' ' ) and t = (ti, ^2, ■ ' ' ) of time 
variables as 



. _Pk_ r _ _Pk_ 

' ~ A; ' ' ~ k 



and consider the generating function 



AGP 



of the new variables. Z[t, t] is a tau function of the Toda hierarchy at a point 
of the lattice P 0- Reversing the sign of t is conventional in the formulation 
of the Toda hierarchy [21] (cf. the fermionic representation of Toda tau 
functions reviewed in Section 3.4). Any Toda tau function thereby becomes 
a tau function of the two-component KP hierarchy. Let us note that the the 
roles of t and i in Z\t,i] can be interchanged by virtue of the identities 

saW = (-l)l^lsu[-t], saH] = (-l)l^lsu[t] (20) 

of the Schur functions and the property ([8]) of ka. 



generalization of this tau function was first studied by Kharchev et al. 
different context. 
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2.4 Cut-and-join operator 

The cut-and-join operator Mq [31] may be thought of as an infinitesimal 
symmetry on the space of tau functions of the KP hierarchy fl2[ [22| [23] . In 
the KP time variables t, the cut-and-join operator reads 

The Schur functions sx[t\ are eigenfunctions of this operator with eigenvalues 
Kx/2: 

MoS,[t] = ^s,[t] (22) 

A combinatorial proof of this fact is presented in Zhou's paper [H]. As we 
shall remark in Section 3, the cut-and-join operator has a fermionic counter- 
part [6], which leads to another proof of ( 122|) . 
(122]) implies the identities 

Therefore one can use e^^'^ to recover Z[t] and from their "initial 

values" at /3 = as 

Z[t] = e^''^Z[t]\p=o, Z[t,t] = e^''^Z[t,th=o. 
Z[t]j3=Q and Z[t,t]|^=o can be calculated by the Cauchy identity 



oo 



J2 sx[t]sx[-t] = exp - ^ ktkh (23) 

AG-P \ k=l / 

and the weighted homogeneity 

sx[cti,cH2,---] =cl^lsA[ti,t2,---] (24) 
of Schur functions as 

Z[t]\^=o = 5^gl'ls,[t]sA[l,0,0, ■ ■ ■] = e<3*i 
xev 

and 

(oo 
k=l 
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One can thus derive the well known formula |12[ [221 ES] 

Z[t] = e^*^Oe«*i (25) 

and its extension 

Z[t, i] = e^"^' exp Q'ktki^ (26) 

to double Hurwitz numbers. 

3 Fermionic representation of tau function 

3.1 Two-dimensional free fermion system 

Let us introduce two-dimensional complex free fermion fields 

ngZ neZ 

Note that we follow the notations of our previous work [321 E3] to use integers 
rather than half- integers for the labels of Fourier modes ■?/'„,'?/'*. The Fourier 
modes satisfy the anti-commutation relations 

The Fock space H of ket vectors and its dual space Ti* of bra vectors are 
decomposed to charge-s sectors Tig, Til, s E Z. Let {s\ and |s) denote the 
ground states in Tig and 7{*, namely, 

{s\ = (-00| ■ ■ ■ V's-lC \s) = 1p-s^^s+l ■ ■ ■ I - oo), 

which satisfy the annihilation conditions 

ipn\s) = for > —s, ipn\^) = for 72 > s + 1, 
{s\ipn = for < — s — 1, {s\ipn = for 77, < s. 

Excited states can be labelled by partitions A = (Ai, A2, ■ ■ ■ , A„, 0, 0, ■ ■ ■ ) of 
arbitrary length as 

|A, S) = ^„Ai-s ■ ■ ■ ijj-X^-s+n-l^t-n+l " " " 

(A, S\ = {S\^_s ■ ■ ■ ^-s+n-li^l^+s-n+l " " " ^Ai+s" 
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|A, s) and (A, s| represent a state in which the semi-infinite subset {Aj + s — 
i + (sometimes referred to as the "Maya diagram") of the set Z of all 

"energy levels" are occupied by particles. These vectors give dual bases of of 
Us and Ti* in the sense that 

{\,r\fi,s) = Sx^,Srs■ (27) 
The normal ordered fermion bilinears 

■.^.irf = ^-ir, - (oi^-,,^;io), i,j e z, 

span the one- dimensional central extension gl(oo) of the Lie algebra gl(c)o) 
of infinite matrices [121112]. gl(oo) consists of infinite matrices A = (ajj)jjgz 
of "finite-band type" , namely, there is a positive integer N (depending on A) 
such that Qij = if \i — j\ > N. For such a matrix A G gl(oo), the fermion 
bilinear 

A = ^ aifip-itfj*: 

becomes a well-defined linear operator on the Fock space, and preserves the 
charge, namely, 

(A,r|l|/i,s) = ifr^s. (28) 

Moreover, for two such matrices A,Bg gl(oo), the associated operators A, B 
satisfy the commutation relation 

[A,B] = \AB]+j{A,B) (29) 
with the c-number cocycle term 

7(A, B) = Tr(A+_i?_+ - i?+-A_+), (30) 
where A±^ and -B±=p denote the following quarter blocks of A, B: 

A^ — = {O'ij)i>0,j<0, = {(lij)i<0,j>0, 

B+- = (^ij)i>0, j<0) = (&ij,)i<o,jr>0- 

3.2 Special fermion bilinears 

The following fermion bilinears are building blocks of our Toda tau function: 

= ^ ■■Ip-n+mlpn-^ m G Z, 
neZ 

Lq = ^n:ilj_nipn-, Wo = ^n^:ilJ-n'ipn-- 

ngZ neZ 
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Jm's span a U(l) current algebra. Lq and Wq are zero-modes of Virasoro and 
W^^^ algebras. These fermion bilinears are associated with infinite matrices 
as 

Jm = ^, Lo = A, Wo = A^, (31) 
where A and A are infinite matrices of the form 

A = (iSij), A = (Si+ij). 
Let J±[t\, t = (ti,t2, ■ ■ ■), denote the special linear combinations 

CXI oo 
k=l k=l 

of Jm's. Their exponentials act on the ground states (s|, |s) as 

(s|e^+W = J2{^,s\s,[t], e^-W|s) = ^ s,[t]|A, s), (32) 
Ae-P \eP 

yielding Schur functions as matrix elements [181 HH] : 

sx[t] = (A,s|e-^+[*l|s) = (s|e'^-W|A,s). (33) 

Unlike other J^'s, Jq is diagonal with respect to |A, s)'s: 

(A, s|Jo|/i,s) = 6xf,s. (34) 

Lq and Wq, too, are diagonal. The diagonal matrix elements can be calculated 
as follows. 



Lemma 1. 



(A, s\Lo\^i, s) = (^|A| + , (35) 

(A, s\Wo\fi, s) = 6,, + {2s + 1)|A| + f(i±lK^i±l) . (36) 
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Proof. Assuming that s > 0, one can calculate the diagonal matrix elements 



as 



(A, s|Lo|A, s) = ^^(Aj + s — i + 1) — ^^(~^ + 1) (heuristic expression) 

1=1 i=l 
oo s 

= ^^((Aj + s — i + 1) — (s — i + 1)) + k (re-summed) 

i=l k=0 
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and 

oo oo 

(A, s| Wol-^; = ^^(Aj + s — i + — ^^(""^ + 1)^ (heuristic expression) 



i=l i=l 

oo 



^(Ai + s-i + lf - {s-i + 1)^) + ^k'^ (re-summed) 

i=l k=0 

s(s + l)(2s + 1) 



Kx + {2s + 1)\X\ + 
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In the case where s < 0, we have only to replace the intermediate sums over 
k as 

k=0 k=s+l k=0 k=s+l 

ending up with the same final expression of the matrix elements. □ 

3.3 Toda Tau function for double Hurwitz numbers 

A general tau function of the Toda hierarchy has the fermionic expression 

T{s,t,t) = (s|e-^+[*l(7e-^-[*l|s), 

where g is an element of GL(oo) [42j. Inserting the aforementioned expansion 
(132|) . one can expand this function as 

Tis,t,t)= ^ {X,s\g\n,s)sx[t]Sf,[-t]. (37) 

Following Okounkov [6|, we now consider the special case where 

g = e^^o/2Q^o = Q^"e^^«/l ^gg) 

Theorem 1. The tau function determined by (E^j can be expanded as 

\ev 

Proof. The properties fl35|) and fl36l) of Lq and Wq imply that g is also diag- 
onalized on the basis |A, s) of the Fock space. The diagonal matrix elements 
take such a form as 

(A, s\g\X,s)= exp (f + (2. + 1)|A| + £(£±IK!£±i) ) ) 

The tau function in question can be thereby expanded as fl39|l shows. □ 



14 



This is a restatement of Okounkov's result [B]. One can rewrite ( 15^ as 

T{s,t,t) = e''^(^+i)(2^+i)/^2Q^(^+i)/2Z^,,.(.+v2,Q[t,t], (40) 

where Z^ q[£, t] denotes the generating function f lTI?]) with the parameters 
f3 and Q being exphcitly indicated. Thus, apart from some numerical fac- 
tors depending on s, the tau function coincides with the generating function 
of double Hurwitz numbers. Note that the s-dependence shows up in the 
generating function as the multiplier e^^'^'^^^'^'^ of the parameter Q. 

Let us conclude this section with a few comments on the cut-and-join 
operator f l2ip . The cut-and-join operator corresponds to the fermion bilinear 

ngZ ^ ^ 

One can readily see from (15^ and (15^ that this operator acts on |A, s)'s 
as 

Mo|A, ^) = (y + ^|A| + |A, (42) 

In particular, |A) = |A,0) is an eigenstate with eigenvalue kx/2. Note that 
s-dependent extra terms show up in the charge-s sector. 

Fermion bilinears of this type can be converted to differential (or "bosonic" ) 
operators by the boson-fermion correspondence [18l[19]. In a generating func- 
tional form, the normal-ordered product 

■.iIj(z)iIj*(w): = ip(z)i!*(w) ^ — (l^;! < Iwl) 

z — w 

of the fermion fields corresponds to the two-variable vertex operator 
X{z,w) 

as 

(s|e^+[*l:^/^(z)^/'*(w): = X(z, w)(s|e^+W. (43) 

A similar relation holds for e'^^'^^'ls) and leads to bosonization with respect 
to i [20], though we omit details here. 
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T{z, w) can be expanded in powers of z — w, and the coefficients of this 
expansion give a bosonic representation of fermion bihnears. For Lq, Jq and 
Wo, this bosonic representation reads 

^o = Y.kt,^+'-^^, Jo = s (44) 

k=l 

and 



k,l=l 



dtkdti 

d s(s + l)(2s + l) 

fc=l 

H) is thus bosonized as 

92 



1 / d 
^o = 7^y2\ kltktiTT^ — + {k + l)tk+l 

2 \ dtk+i 



dtkdti 



Ed As^ — s , , 

fc=l " 



In the charge-0 sector, this reduces to the cut-and-join operator (pTi) . 

4 Generalized string equations for double Hur- 
witz numbers 

4.1 Notations for difference operators 

Building blocks of the Lax formalism of the Toda hierarchy are one-dimensional 
difference operators in the lattice coordinate s [21] . Those operators are lin- 
ear combinations of the shift operators e"^% e"'^^f{s) = f{s + n). Although 
a genuine difference operator is a finite linear combination 

N 

A = ^ a„(s)e"'^^ (operator of [M, iV]-type) 

n=M 

of the shift operators, one can consider a semi-infinite linear combination of 
the form 

TV 

A= ^ a„(s)e"^" (operator of (-00, A^] type) 
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or 

oo 

A = ^ a„(s)e"^'' (operator of [M, oo) type) 

n=M 

as well, which amount to pseudo-differential operators in the Lax formalism of 
the KP hierarchy [18]. Let ( )>o and ( )<o denote the truncation operation 

(A)>o = 5^a„(s)e"^% (A)<o = 5^a„(s)e"^^ 

n>0 n<0 

Difference operators are in one-to-one correspondence with ZxZ matrices. 
Firstly, the n-th shift operator e"^'' corresponds to the shift matrix 

= {^i,j-n)i,jeZ- 

Secondly, the multiplication operator a{s) amounts to the diagonal matrix 

diag(a(s)) = (a(05ii)i,jGZ- 
In particular, the multiphcation operator s corresponds to 

A = diag(s) = ii5ij)ijez- 
Consequently, a general difference operator 

n 

is converted to the infinite matrix 

A{A,A) = ^diag(a„(s))A" = J^(an(«)'^i,i-n)i,iez- 

n n 

Occasionally, it might be more convenient to write a difference operator in 
an anti-normal-ordered form as 

i?(e^%3) = 5^e"^^6„(s). 

n 

In that case, the corresponding infinite matrix reads 

S(A,A) = ^A"diag(6„(s)). 

n 
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4.2 Lax and Orlov-Schulman operators 

The Lax formalism of the Toda hierarchy uses two Lax operators L, L of 
type (—00, 1] and [1, 00). Actually, from the point of view of symmetry, it is 
better to consider L and rather than L and L. These operators admit 
freedom of gauge transformations. In the gauge where L is monic (namely, 
the leading coefficients is equal to 1), L and can be expressed as 

00 

L = e^^ + ^M„e(^-")^% 

n=l 

00 
n=l 

The coefficients Un and till functions of s and the time variables t, t, and 
written as Un{s,t,i) and Un{s,t,i) if we do not suppress the independent 
variables. L and L satisfy the Lax equations 

— = [Bn,L], -^ = [i?„,L], 

-\B n -\B n 

where Br, and are defined as 



Bn = {L^)>Q, Bn = {L "'] 



<o- 



To formulate the generalized string equations, we need another pair of dif- 
ference operators M,M, namely, the Orlov-Schulman operators [H]. These 
operators, too, satisfy the Lax equations 

dM dM 

dM _ dM ^ ' 



[i?„,M], -^ = [Bn,M] 



dtn ' ' dtn 

of the same form as the Lax operators do, and are related to the Lax operators 
by the twisted canonical commutation relations 

[L, M] = L, [L, M] = L. (49) 

By "generalized string equations" we mean equations of the form 

C{L,M) = C{L,M), (50) 

where C{L,M) and C{L,M) are (possibly infinite) linear combinations of 
monomials of L, M and L, M with constant coefficients. The following lemma 
[30] l3T] explains an origin of generalized string equations. 
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Lemma 2. If the fermion hilinears C(A, A) an(iC(A, A) are intertwined by 
a GL(oo) element g as 

cXK^)g = gClj^), (51) 
then the Lax and Orlov-Schulman operators satisfy /f5Q|). 

4.3 Intertwining relations 

Let us now consider the case of the solution determined by the GL{oo) ele- 
ment f l38|) . We seek intertwining relations in the following special form: 



Jkg = gC{A,A), C{A, A)g = gJ.k, k = l,2,---. 
Lemma 3. Jm 's are transformed by the adjoint action of and e^^^l"^ as 

Q JmQ Q Jmi 

jiez 

Proof. Let us note the fundamental commutation relations 
[Lo, ilJn] = -mpn, [Lo, ^A*] = -mjj^ 

and 

that follow from the definition of Lq and Wq- These commutation relations 
can be exponentiated as 

and 

Consequently, the exponentiated operators Q^" and e^^°^'^ act on the basis 
■.ip-iip*: of gl(oo) as 

and 
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One can thereby derive (1^ as 



nGZ 

and 



neZ 



neZ 

□ 

Lemma 4. J±fc 's are connected with the fermion bilinears A^e^^'^ and X-^e^^^ 
by the GL(oo) element ( Tfg)) as 

Jfc^? = gQ^e-^^''^K^^^, (53) 
^J_fc = Q'^e^^^'^I^^^^g. (54) 

Proof. Using the relations ( 152|) in the previous lemma, one can calculate 
5f"Vfe5f as 

neZ 

Since the last sum can be rewritten as 

nGZ 

the first intertwining relation ( 153|) follows. In the same way, one can calculate 
gJ-kg~^ as 



nGZ 

which implies the second intertwining relation f l54p . □ 
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4.4 Generalized string equations 

Theorem 2. The Lax and Orlov- Schulman operators of the tau function 
satisfy the generalized string equations 

= Q^e-^'^'/^L'^e^''^', = Q^e^^" /^L-^e^^^ (55) 

for k = 1,2, ■ ■ ■ . Moreover, these equations can be derived from the first two 
(k = 1) equations 

L = Qe-^/^Le^^, L~' = Qe^/^L-'ef'' . (56) 

Proof. The first part is a consequence of (1531) and (I54p . Let us show the 
second part. The fc-th power of the first equation of fl56|l reads 

The commutation equation of L and M in implies that 

[M, ■ ■ ■ , [M, [M, I]] ■ ■ ■ ] (fc-fold commutator) = {-Vfh 
for /c = 1, 2, ■ ■ ■ , so that 

^mi^-m = L + V ^[M, ■ ■ ■ , [M, [M, L]] • ■ ■ ] = e-f^L. 
k=i ^■ 

Using this relation repeatedly, one can move e'^^'s in 

(^^g/3M)fe ^Yie right- 
most position as 

^ g-/3fc-/3(fc-l)_^3g/3M^g/3Mj^^A:-3g3/3M 

Thus the first equation of (l55i) follows. The second equation of (l55|) . too, can 
be derived from (156|) in the same way. □ 

Thus, in contrast with two-dimensional quantum gravity [25] and c = 1 
string theory [2Sl [23 [2S1 123 [SDl 131], the generalized string equations contain 
the exponential terms e^^ , e^^ . These terms stem from the fermion bilinears 
yY±fcg/3fcA f[53|) and fl54|) . Ferimion blinears of a similar form are also used 
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in the study of integrable structures of the melting crystal model [221 E3]- 
A common algebraic background of these fermion bilinears is the quantum 
torus algebra (with parameter q) spanned by the infinite matrices 

= q-km/2^m^k/^^ k.TJieZ. (57) 

In the melting crystal model [521 [33], a central extension of this Lie algebra 
is realized by the fermion bilinears 

Ki'^ = (58) 

neZ 

The Lax and Orlov-Schulman operators give (two copies of) yet another kind 
of realization by the difference operators 

ylj.) ^ q-km/2j^m^kM (gg) 

and 

y(k) ^ ^-km/2lm^kM_ (gg) 

If q = e^, the exponentials e^*^, e^*^ belong to this Lie algebra. 

5 Classical limit of generalized string equa- 
tions 

5.1 Dispersionless Toda hierarchy 

In a naive sense [43], the classical limit of the Toda hierarchy can be ob- 
tained by replacing the shift operator e*^* by a new variable p. The difference 
operators L, M, L, M thus turn into Laurent series of p of the form 

oo 

C = p + ^UnP^'"", 
n=l 

oo 

£-1 =Uop-^ + J2unP''~\ 
n=l 

oo oo 

n=l n=l 

oo oo 

n=l 71=1 
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that are referred to as the "Lax and Orlov-Schulman functions" . 

As difference operators arc replaced by Laurent series, commutators of 
difference operators turn into Poisson brackets by the rule 

[e^%s] = e^' {p,s} = s. 

Accordingly, Poisson brackets of functions of p and s are defined as 



{F,G}=pi 



^d£dG _ d£dG 
\ dp ds ds dp 



The Lax equations and the twisted canonical commutation relations are re- 
defined with respect to the Poisson bracket as 

djC dC 



dM dM 

Otn Ot 



(61) 



n 



and 

{£,M} = A {C,M}=C. (62) 
Bn and Bn are defined by seemingly the same formulae 

Bn = (>C")>o, Bn = {C ")<o 

as in the previous case, but the notations ( )>o and ( )<o now stand for 
projection operators on the space of Laurent series, namely, 

X]«nP"J =^anP'', \^(lnf"\ =^anP''. 
n / >0 "^^0 \ " / <0 "<0 

These equations are fundamental constituents of the "dispersionless Toda 
hierarchy" . 
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5.2 /i-dependent Toda hierarchy and generaUzed string 
equations 

The foregoing procedure replacing e^" — )• p can be justified as a kind of 
classical limit in an /i-dependent formulation of the Toda hierarchy [H] . 

In the ^-dependent formulation, e^" is replaced by e''^" . The "Planck con- 
stant" h thus plays the role of lattice spacing. The Lax and Orlov-Schulman 
operators are expanded in powers of e^^" as 



-n)hdB 
n=l 



oo 
n=l 

oo 



n=l 

oo 
n=l 

oo 

M = ^ nt„L" + s + ^ Vnl-"", 

n=l n=l 
oo oo 

M = - ^ nt„L-" + s + ^ Vnl"". 

n=l n=l 

The Lax equations and the twisted canonical commutation relations take an 
fi-dependent form as 

h— 

dtn 

dl 



^dtn 
dtn 

dM 

dtn 



h 



■ [Bn,L], 


h— = 

dtn 


[Bn-, L]-, 


'- [Bn,L], 




[Bnj L], 


[Bn,M], 


din 


^ [S„,M], 




^dM 
dtn 


= [Bn,M] 



(63) 



and 



[L, M] = hL, [L, M] = hL. (64) 



If the coefficient m„, Un-, Vn, Vn (which are functions of h, s, t, t) have a 
smooth classical limit as 

u^^ = limti„, 'ui°'' = lim-Un, v^^ = limt>„, vi^^ = lim-u^, (65) 
h^o Ti-i-o h^o h^o 
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one can define the Lax and Orlov-Schulman functions as 

oo 
n=l 

oo 
n=l 

oo oo 

n=l n=l 

oo oo 

M = -Y^ ntnC.-'' + s + Y^ 

n=l n=l 

Tliese Lax and Orlov-Scliulman functions satisfy the Lax equations (16T|) and 
the twisted canonical Poisson relations (l62l) . 

In this /i-dependent formulation, generahzed string equations (150|) are 
modified as 

/i^^M) = C(L, n^^M), (66) 

namely, M and M are multiplied by [201 EI]- Let us explain the underly- 
ing mechanism briefly. The Lax and Orlov-Schulman operators are connected 
with the matrices A and A by the so called "dressing operators" [44]- The 
twisted canonical commutation relations are thereby derived from the com- 
mutation relation 

[A,A] = A 

of these matrices. In the /i-dependent formulation, they take the form 

[L, h-^M] = h-^M, [L, h-^M] = h-^M, 

which are nothing but flMl) . Thus it is h~^M and h~^M rather than M and 
M that correspond to A and show up in generalized string equations. 

5.3 Classical limit of generalized string equations 

Let us turn to the case of double Hurwitz numbers. To derive a classical 
limit, we have to introduce h therein. At least formally, this can be done by 
rescaling the space-time variables as 

s h-^s, t h'H, i frH. (67) 
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The /^-independent Toda hierarchy is thereby converted to the /i-dependent 
form. Actually, it is rare that the rescaled Lax and Orlov-Schulman operators 
satisfy the condition (1651) . To achieve a meaningful (and nontrivial) classical 
limit, one should start from a carefully chosen h- dependent solution of the 
h-independent Toda hierarchy. 

In the language of the tau function [44], a meaningful classical limit is 
obtained from an /i-dependent tau function T{h, s, t, i) such that the rescaled 
tau function 

Tt{s, t, i) = T{h, h~^s, h'H, h'H) 

behaves as 

\ogTn{s,t,t) = h-^J^{s,t,i) + 0{h-^). (68) 

J-" = J^{s, t, i) is called the "free energy" because of its relation to random 
matrices and topological field theories [IS, |16] . If the rescaled tau function 
has this asymptotic form, the associated "wave functions" 

»>.(..M,^)= "'''''-'r':'^:''-^' .-'"-e'-^"'". 

Tn{s,t,t) 

This, t, t) 



k=l 

have the "WKB" form 

^r,,(s, t, t, z) = exp {h-^Sis, t, t, z) + 0(1)) , 
^!f,{s, t, i, z) = exp {h~^S{s, t, i, z) + 0(1)) 



(69) 



and satisfy a set of auxiliary linear equations. The phase functions S'(s, t, t, 2;) 
and S{s,t,i,z) satisfy the associated Hamilton- Jacobi equations, which can 
be converted to the dispersionless Lax equations (!6T|) and the Poisson rela- 
tions (see the review [20] for details). 

An appropriate /i-dependent reformulation of the tau function (15^ of 
double Hurwitz numbers can be found by the following heuristic considera- 
tion. As we move into the /i-dependent formulation, the generalized string 
equations (l56l) are modified as 
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Obviously, these equations do not have a hmit as /i — )■ 0. If, however, the 
parameter (3 is simultaneously rescaled as 

/3 h(3, (70) 

the generalized string equations are further modified as 

L = Qe-^^/2Le'3^^ = Qe^^/^^-^e^^, (71) 

and have a meaningful classical limit of the form 

C = QCe^^, C-^ = QC-'e^^. (72) 

5.4 Existence of /i-expansion 

To justify the foregoing heuristic derivation of the classical limit (1721) of the 
generalized string equations, let us show that the tau function (15^ with /3 
rescaled as fITOl) does satisfy the condition fl68l) . 

Recall the expression fHOj) of the tau function. After rescaling s, t, i and 
/3 as fl671) and flTOj) . this expression is modified as 



nis, t, i) = e^-'''^(^+'')(2^+'^)/^2g'"'"'^(^+^)/2Z,^,,.(..V2)Qr 4, h-% (73) 

Therefore it is sufficient to show that the logarithm of Zfif^^qlh^^t, h^H] has 
an ^-expansion of the form 

log Znp,Q[h-H, h-H] = h-^Fo + Fi + h^F2 + --- + te-^Fn + • • • , (74) 

where Fq, Fi, ■ ■ ■ are analytic functions of (/3, Q, i) in a common domain. 
The free energy is then given by 

•F=^ + ^^ + W,e^^Q,t,t). (75) 
o 2 

(I71|) is a generalization of the well known topological expansion for simple 
Hurwitz numbers. It is common in the literature that this kind of expansion 
is explained by a combination of topological and combinatorial consideration 
(see, e.g.. Section 4.2 of Mironov and Morozov [22j, Section 2.1 of Bouchard 
and Marino [13] and Section 2.2 of Borot et al. [15J). We take another 
approach based on the cut-and-join operator fl?T]) . 

Theorem 3. log Zrii3^Q[h'^t,h^^t\ has an h-expansion of the form [74^ . 
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Proof. Let F = F{h, (3, Q, t, t) denote the left hand side of ( 1711) multiplied 
by /j.^. By fl2B]) . e'^ ^ can be expressed as 

(CXI 
k=i / 

where Mo{h) denotes the rescaled cut-and-join operator 

92 



1 / d 

2 V otk+i 



J, 



dtkdti 



Therefore e satisfies the differential equation 



hMQ{h)e- 



dp 

with respect to /3. This equation can be further converted to the differential 
equation 

d'^F dF dF 



dF 1 dF 1 ( 



+ 



9tfc(9ti dti 

for F. This equation is supplemented by the the initial condition 

00 

-F|/3=o = ~ Q'^ktkik- 

k=l 

We now seek a solution of this initial value problem in the form of a (formal) 
power series of h'^: 

F = Fo + h^Fi + --- + h^'^Fn + ■■■ . 
This reduces to solving the differential equations 

dF^ 



' ki=i ki=i ' 



kl=l m=0 ' 



for n = 0, 1, 2, ■ ■ ■ under the initial conditions 



n|/3=0 = 

k=l 



5noX]Q'A;44. (77) 
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FnS are thereby recursively determined, and become analytic functions of 
{P,Q,t,i) in a common domain of definition (because the differential equa- 
tions other than the first one for n = are linear with respect to Fn). Since 
the initial value problem for F has a unique solution, the power series solution 
F = Fo + h^Fi H should coincide with the left hand side of □ 



One can thus confirm the expected asymptotic form fl68l) of the the 
rescaled tau function fl73|) . Let us stress that this is also a proof for the 
case of simple Hurwitz numbers. To consider that case, one has only to set 
ik = Ski in the initial condition fl77|) . Let us also point out that the main 
part Fq of the free energy is determined by the n = part of ([76]) and (177j) . 
namely, 

k,l=l k,l=l " ' 

and 

Fo|;3=0 = -5^Q''Hfctfe. (79) 
k=l 

It will be interesting to apply the diagramatic technique of Mironov and 
Morozov [22] to these equations. 



6 Solution of classical limit of generalized string 
equations 

6.1 Comparison with generalized string equations of 
c = 1 string theory 

Our goal in this section is to solve the generalized string equations ( 172]) and 
to derive some implications thereof. To this end, it is instructive to compare 
these equations with the generalized string equations of c = 1 string theory 

[2Sl 127] [221123. 

In the classical limit, the generalized string equations of c = 1 string 
theory read 

C = CM, C-^ = C-^M. (80) 

Let us mention that the same equations play a central role in a problem 
of complex analysis and its applications to interface dynamics |17l [H] 
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Actually, it is the equation 



= 1 (81) 

rather than flHOj) that is referred to as a "string equation" in these applica- 
tions. Note that one can readily derive flHTj) from flHOj) . In the same sense, 
one can derive the equation 

{log£,log£-i} = /3 (82) 

from fl72|) as a counterpart of flHOl) for double Hurwitz numbers. 

flHTj) and fl82|) resemble the string equation (or the Douglas equation) 

[Q,P] = 1 (83) 

and its classical limit 

{g,p} = i (84) 

in two-dimensional quantum gravity [501 EB |52] • Q and P in (183|) are one- 
dimensional differential operators of the form 

Q = dl + Ma^r^ + ■■■ + «„, P = d': + V2d^-^ + --- + Vn. 

In the classical limit, dx is replaced by a variable p with the Poisson bracket 

= 1, 

and Q and P are polynomials of the form 

Q = + U2P"-' + ■ ■ ■ + Un, P = + V2P'^-^ + --- + Vn. 

In this setting, Q and P may be thought of as coordinates of the spectral 
curve (parametrized by p) in the sense of Eynard and Orantin [l7j. Namely, 
when X and other deformation variables (time variables of the underlying KP 
hierarchy) are fixed to special values, Q and V satisfy a defining equation 
/(X, y) = of the spectral curve as 

/(Q,^) = o. 

Although the KP and Toda hierarchies are different in nature, the last 
remark seems to suggest that one may think of an equation of the form 

/(£,£-i) = 
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as the spectral curve in the present setting. This observation is partly sup- 
ported by the fact that such a curve is derived as the spectral curve in 
the random matrix approach to c = 1 string theory and interface dynamics 

[531 EH Eg. 

Bearing these remarks in mind, let us turn to the issue of solving the 
generalized string equations ( 172]) . These equations, like ( 15U]1 . are a kind of 
"nonlinear Riemann-Hilbert problems". One can use the method developed 
for solving flHOj) [29] to construct a solution of fl72l) as power series of t and 
i. As it turns out, fl72|) can be treated in much the same way apart from 
technical complications. 

6.2 Decomposition of equations 

Let us convert fl72|) to the logarithmic form 

log(£p-i) = logg - log(£-V) + (3M, 



log(£-V) = logg - \ogiCp-') + (3M. 
Since Cp~^ and C^^p are Laurent series of the form 

oo oo 
n=l 71=1 

with nonzero leading terms, one can expand the logarithm as 

oo oo 

log(£p"^) = a„p"" , log(£"V) = log Mo + ^ ttnP" 

n=l Ti=l 

where 

c^n = Un + (polynomial of Mi, ■ ■ ■ , m„_i), 
ttn = UQ^Un + (polynomial of Uq^ui, ■ • • , UQ^Un-i)- 

We now substitute 

oo oo 

k=l n=l 

oo oo 



k=l n=l 



15) 



in fl85l) and expand both hand sides in powers of p. This leads to an infinite 
set of equations for the coefficients Un, Un, fn, Vn of C, Ai, C, M. as follows. 
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Equating the coefficients of p = 0, 1, ■ ■ ■ , in both hand sides of the 
ffist equation of f lS^ gives the equations 

CXI 

= log g - log Mo + /3s - /3 5^ kik{C'%, (86) 

fc=i 

oo 

an = -PY,kik[C-'')-n, n = l,2,---, (87) 

k=n 

where {C~^)^n stands for the coefficient of in C^^ . In the same way, 
the coefficients of p", = 0, 1, 2, ■ ■ ■ , in the second equation of fl85|l give the 
equations 

cx> 

\ogu^ = \ogQ + l3s + pY.ktk{C%, (88) 

k=l 

oo 

a„ = /3^Hfc (£'=)„, n = l,2,---, (89) 

k=n 

where {C^)n denotes the coefficient of in C^. Since 

{C% = {C-')-n = for A; < n, 

the range of k in the sums of (|57|) and is limited to /c > n. Note that 
one can use the formal residue notation 




= a_i 



to express {'C.'')^ and (£ '^)_^ as 

{C% = res (£V"^^logp) , (£-'^)-„ = res (£-Vt/logp) . 

Such an expression turns out to be useful in the subsequent consideration. 
Since VnS and VnS are absent, ( 15^ - (jHSD are equations for m„'s and Mb's 
only. 

The remaining part of fl85|l determine f„'s and w„'s. To see this, it is more 
convenient to expand fl85|) in powers of C and £ rather than of p. Extracting 
the coefficients of from the ffist equation of fl721) and those of from 
the second equation yields the equations 

= - res (log(£- V)>C""c/ log C) + (3vn, (90) 
= - res (log(£p~^)£"ci log C) + (3vn (91) 
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for n = 1, 2, ■ ■ ■ . Thus f„'s and f„'s are determined by £ and £ as 



Vn = /3-^res (log(£p-^)£"rflog£) , 
Vn = p-^res {log{C-^p)C-''dlogC) . 



( 185|) can be thus decomposed into the infinite set of equations (l86i) - (189|1 . 
(190|) and (l9Ti) . The next task is to show that they do have a solution. 

6.3 Solution of equations 

As we observed above, one can think of fl86l) - fl89l) as equations for m„'s and 
UnS. We want to construct these functions as power series of {t and i with 
coefficients depending on s. If UkS and -Ufc are thus constructed, Vk and v^s 
are determined by fpJ]) and f PT]) . 

Let us first examine (ISBj) and f lSS]) . Since 

(yn = Un + (monomials of higher degrees in mi, • • • , ti„_i), 
«n = Uo^Un + (monomials of higher degrees in Uq^u^ ■ ■ ■ , Mq ^-Un-i), 

Un and u„, = 1, 2, ■ ■ ■ , show up on the left hand side of these equations 
linearly. The right hand side consists of terms that are multiplied by t^'s and 
tfc's. Therefore these equations yield a huge system of recursion relations for 
the coefficients of power series expansion of m„'s and UnS. 

Note that uq, which remains to be determined, is contained in the coeffi- 
cients of the power series expansion of m„ and «„. We need another equation 
to determine Uq. Actually, there are two equations and (1551) rather than 
just one. This puzzle is resolved as follows 0. 

Lemma 5. // [8l\ ) and are satisfied, then ^8^) and are equivalent, 
and reduces to the equation 



Proof. Let us use the formal residue notation to express the terms {C^)o in 



similar result holds for the generalized string equations ((80)) of c = 1 string theory. 
This fills a logical gap left in our previous paper [29] . 



oo 




(92) 



k=l 



(186]) as 



(£'^)o = res(£*^(ilogp). 
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Since the identity 



= res(£'^(ilog£) = res ( £V — d\ogp 



dp 

holds for /c > 1, this expression of {CJ')q can be further rewritten as 
(£^)o = res (^C\l ~ p^^^)d\ogp 

,^aiog(£p-i) 



= — res £ dp 

\ op 

By substituting Cp^^ = aip^^ + a;2P~^ + " " " ) the right hand side can be 
expanded as 

RHS = ai ms{C^p~^dp) + 2^2 res(£ V^^^p) + ■ ■ ■ 

Since {C^)n = for n > k, this expansion terminates at the k-th term. One 
can thus obtain the identity 

iC% = a,{C% + 2a2{C% + ■■■ + kak{C%. 

In the same way, one can derive the identity 

{C-% = + 2a2(^-')-2 + • • • + kak{C-^)^k. 

By virtue of these identities, one can rewrite the two sums in fl86l) and 



as 

oo 



k=l k=l 

oo 



^ ktk{C% = ktk + 2a2{C% + ■■■ + kau{C%) 

k=l 

oo oo 

«i ktk{C% + 2a2 Y k^k{C% + ■■■ 

k=l 

oo 



fc=i fc=i 

oo 



n=l 

and 

oo oo 

Y kh{C'% = Y kh {o^i{C-% + 2a2{C-% + ■■■ + kak{C-%) 

k=l k=l 

oo oo 

= «i 5^ kh{C-% + 2^2 Y kh{C'% + ■■■ 

k=l k=l 

oo 
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n=l 



Note that fl5^ and (I57|) have been used to derive the last hnes. Thus 
and fISS]) turn out to reduce to the same equation □ 

We can thus use fl^ in place of fISB]) and f lSS]) . Adding this equation to 
(187|) and (l89|) . we obtain a full system of equations that determine the power 
series expansion of Un and uo recursively. We can readily see from ( 192|) 
that Mo is a power series of the form 

log Mo = logQ + /3s + (terms of positive orders in t, i). (93) 

On the other hand, since 

(£")„ = !, {C-n-n = ul 
{C^)n = (polynomial in -Ui, ■ ■ ■ , Uk-n) for k > n, 
{C^^)^n = ^ (polynomial in Mq^Mi, ■ ■ ■ ,UQ^Uk-n) for k > n, 

Un and Un are power series of the form 

Un = —f^ninUQ + (terms of higher orders in t, i), ^^^^ 
Un = f3ntnUo + (terms of higher orders in t, i). 

This power series solution of fl72|) (which is unique by construction) is 
homogeneous just like the solution of the generalized string equations flHOjl 
for c = 1 string theory |29]. This is a consequence of invariance of the the 
string equations under the scaling transformations 

tn C'^tn, in cHn, S ^ S, p ^ Cp, ^^^^ 
Un c"m„, Un C""m„, Vn C^Vn, Vn C""t;„. 

In summary, we have observed the following: 

Theorem 4. The generalized string equations ( [7^ have a unique solution 
that has power series expansion with respect to {t, i) as shown in and 
[94\ )- This solution is homogeneous with respect to the scaling transformation 

6.4 Solutions at special values of t, i 

The foregoing construction of solution simplifies to some extent when t and 
i take special values. Of particular interest are the following two cases: 

(i) tfe's are free, and ikS are restricted to 4 = ii6ki, 

(ii) tfc's are restricted to ik = ii^ki, and t^s are free. 
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They amount to restricting the generating function Z[t, i] of double Hurwitz 
numbers to generating functions of simple Hurwitz numbers. Since these two 
cases are essentially equivalent, let us consider (i) only. 

In the case of (i), (187|) implies that an vanishes for n > 1 and that the 
only non-vanishing component is given by 

«! = Ul = —PiiUQ. 

Thus £ simphfies as 

C = pe'^'P'' = pe-^*^"0P"' , (96) 

and (£'^)„ can be written explicitly as 

, (fc^i)^-" (-fc/3ti^o)^- 
^ ^" ~ ik-n)\ ~ ik-n)\ ■ ^ ^ 

Un, n = 1, 2, ■ ■ ■ , are thereby recursively determined by (!89|) as a function of 
tfc's and Uq. Uq is determined by (1881) . which now takes an explicit form as 

log^o = log Q + /3. + /3 kt.thl^p^. (98) 

k=l 

Vn and VnS, too, have more or less explicit formulae, though we omit details. 

This result shows a remarkable feature. Namely, ( 196|) resembles the defin- 
ing equation x = ye^ of Lambert's W-function y = W{x) if and are 
identified with x and y. The so called Lambert curve is defined by this equa- 
tion on the (x, ?/)-plane, and plays a fundamental role in the recent studies 
on Hurwitz numbers [HI HSl [ISl ESI ES| • 

This analogy becomes more precise when t^s, too, are specialized to 
tk = 0, k = 1,2, ■ ■ ■ . In that case, uq is explicitly determined as 

uo = Qe^'. 

Moreover, (l89l) implies that an vanishes for all n, hence 

£-1 = uop-\ (99) 
fl96|l thereby turns into the equation 

C = uoCe-^''^~' (100) 

for C and £. In view of the remarks in the beginning of this section, it seems 
likely that this equation can be identified with the spectral curve for simple 
Hurwitz numbers. 
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